Low-temperature resistivity saturation in the Kondo insulator SmB 6 : 
Indication for another scenario of hopping-type conductivity 
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We provide arguments that dynamical changes of energies of impurity states due to valence 
fluctuations (VFs) result in non-equilibrium energy distribution function of localized charge carri- 
ers. Corresponding hopping type of conductivity that does not require a thermal activation can 
explain unusual metallic-like conduction of SmBg and other Kondo insulators experimentally ob- 
served at lowest temperatures. A qualitatively similar hopping process is expected in metal-oxide- 
semiconductor (MOS) structures due to local dynamical changes of voltage potential in the metal 
(gate) electrode. 

PACS numbers: 75.30.Mb, 72.20.Ee, 72.15.Rn, 73.40.Qv 



SmBg belongs to the class of intermediate valence (or 
homogenously mixed valence) semiconductors™ which 
are also known as Kondo insulators or heavy-fermion 
semiconductors^. The ground state of Sm is coherent 
superposition of the configurations, 4/ 6 and 4/ 5 5d, with 
weight factors 0.3 and 0.7, respectively, giving an effec- 
tive 4/ valence of 2.74 Although this compound has been 
studied for more than four decades, many fundamental 
aspects, especially those associated with the origin of the 
gap and the unusual low-temperature transport proper- 
ties remain still unclear. A principal problem represents 
a metallic-like electrical conductivity of SmBg at lowest 
temperatures. 

Electrical resistivity measurements of SmBg show a 
large resistivity increase at decreasing temperature below 
50 K with a saturation at high residual value po at lowest 
temperatures ) » 5 ' 6 ' 7 ^ Detailed studies of electrical con- 
ductivity <j(T) show a temperature non-activated (metal- 
lic) channel in <r(T), which becomes dominating below 
about 3 K4&1&2 An explanation of the high value of po 
requires a superunitary scattering^ with unphysically 
high concentration of scattering centers (at least 80 per 
unit cell 6 ). However, according to the Mott-Ioffe-Regel 
viewpoint, the conventional Boltzmann transport theory 
becomes meaningless when the characteristic mean free 
path of the itinerant conduction electrons becomes com- 
parable to, or less than, the interatomic spacingj 10 ' 11 ! 12 
Therefore, the high value of the residual resistivity can- 
not be attributed to metallic conductivity mediated by 
itinerant electrons. On the other hand, assuming that 
electron states in the vicinity of the Fermi level are lo- 
calized, the present theories require a temperature acti- 
vated conductivity and an insulating ground state, what 
is, however, in qualitative disagreement with experimen- 
tal indicationsiii£i£i£i£ mentioned above. The consequence 
of this disagreement is that the origin of the residual re- 
sistivity/conductivity in SmBg can not be consistently 
explained considering the known mechanisms of electri- 
cal conductivity in metals or semiconductors. 

In this paper we provide arguments that dynamical 
changes of energies of localized (impurity) states, which 
have to be present in SmBg due to valence fluctuations 
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FIG. 1: Schematic depiction of energy of donor levels in the 
IB of classical semiconductor (a) with corresponding DOS di- 
agram (b), time evolution of the energy of donor impurity 
due to rearrangements of metallic ions (c) and time-averaged 
partial DOS (d) corresponding to SIB (see definition in the 
text). 



(VFs), result in a non-equilibrium energy distribution 
function of localized charge carriers. As a consequence, a 
temperature nonactivated hopping (tunneling) between 
localized states occurs yielding metallic-like conductivity 
at lowest temperatures. We also advert to the fact that 
qualitatively analogous mechanism should be present in 
MOS based devices, e.g. in MOS field effect transistors 
(MOSFETs) because of local dynamical changes of volt- 
age potential in the gate (metal electrode). 

Let us start considering a classic semiconductor with 
a donor impurity band (IB) located in the forbidden 
gap. Let the Fermi energy, Ep, lies in the IB and the 
states near Ep are localized, characterized by a constant 
density-of-states (DOS) function, g(E). In addition, let 
us suppose that the crystal lattice of the semiconductor 
contains metal ions Me in two different valence states at 
least, say Me 2+ and Me 3+ , randomly distributed over 
the lattice. The energy level diagram and the DOS of 
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the system in a "static" case, when the metal ions do not 
change their valence, is schematically depicted in Fig. [TJa, 
and Fig. [TJd, respectively. The electrical conductivity at 
lowest temperatures in such a system is of a variable- 
range hopping typ o 12 i 13 , arising from the phonon as- 
sisted hopping in a so-called optimal band. The opti- 
mal band is a narrow energy interval of the width 2eo 13 
centered at the Fermi level with concentration of impuri- 
ties N(e ) — 2e g(Ep)M- Several derivations consistently 
show that eo, being in fact a typical activation energy of 
a hop, decreases with temperature as T 3 / 4 , causing a 
consequent decrease of the impurity concentration in the 
optimal band. As a result, an average distance between 
the hopping centers increases as T -1 / 4 , and the conduc- 
tivity follows Mott's law 12 ' 13 ' 14 , a oc e ( - To/T)1/4 . 

Now, let us consider such rearrangements of metal-ions 
valences (with a characteristic time constant, i r , suffi- 
ciently long for well defined energies of impurity states) 
that do not change neither the number of Me 2+ and 
Me 3+ ions, nor the physical properties of the system. 
Due to different physical properties of Me 2+ and Me 3+ 
ions (e.g. ionic radius or charge), every such rearrange- 
ment causes a change of local parameters influencing 
the energy of impurity states (e.g. chemical pressure 
acting on the impurity center or the Coulomb interac- 
tion between the impurity center and the metallic ions). 
Consequently, the energies of individual donor states are 
changed too, as it is schematically shown in Fig. [IF,. Per- 
forming all the possible rearrangements as defined above, 
the energy of a donor level, Ei, satisfies the inequality 



< Ei < Ei. 



(1) 



where Ei^ m i n and Ei^ max define the energy window for 
the donor level. So that, because of VFs, the energy 
of impurity state at a site i dynamically changes in the 
interval of the width 



AE, = Ei, 



Ei 



(2) 



For purposes of this paper we denote this interval as a 
single impurity band (SIB). The probability of finding an 
impurity state i at the energy E from its SIB is charac- 
terized by partial-DOS function gi(E) coupled with g(E) 
by 



g{E) = Y,9i{E), 



(3) 



where the summation is performed over all impurity 
states. Such as Ep lies in the IB, there have to exist 
also SIBs satisfying the condition 



Ei 



<E F < Ei. 



(4) 



For these states there is a non-zero probability that some 
occupied donor levels from the region below Ep shift due 
to the valence rearrangements to the region above Ep. 
Analogously, empty donor levels from the region above 
Ep can shift under Ep. Such energy changes represent 
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FIG. 2: Schematic depiction of hoppings due to a rearrange- 
ment within impurity subnetwork given by inequality (4). 
The states before rearrangement are represented by dotted 
symbols, the width of SIBs is shown as "error bars". The 
hopping from occupied state i shifted above Ep is to state j 
directly, or via state k, which is temporarily included in the 
subnetwork. Hopping to the empty state m which decreased 
below Ep from occupied state I is another example of possible 
hopping to the state with less energy. 



in fact non-equilibrium excitations driven by the rear- 
rangement process (RP). Therefore, they have to be con- 
sequently brought into the equilibrium state, namely via 
electron hops from occupied states above Ep to empty 
states with less energy, or via hops to empty states below 
Ep from occupied states with higher energy. For simplic- 
ity let us assume that all SIBs are characterized by an 
equal width AEi w Eq . In such case energies Ei from the 
SIBs satisfying the inequality (4) lie in the band given by 



Ep — Eq < Ei < Ep + Eg. 



(5) 



The corresponding donor states, as well as those only 
temporarily satisfying the inequality (5), form a subnet- 
work, in which conditions for hops from occupied states 
to empty states with less energy are intrinsically created, 
as schematically shown in Fig. 2. Because such hopping 
processes do not require activation energy, we denote the 
interval defined by eq. (5) as zero-activation-energy band 
(ZAEB). The number of states in ZAEB is 



N(E ) = 2E g(E F ). 



(6) 



Although Eq.(6) has formally the same form as the above 
mentioned equation for the concentration of hopping cen- 
ters in the optimal band, considered in the variable-range 
hopping conductivity^ 3 -, there is a fundamental differ- 
ence: while the width of the optimal band (2eo) decreases 
with temperature as T 3//4 yielding the zero concentration 
of hopping centers for T — > K, the concentration of 
hopping centers in ZAEB is always non-zero because of a 
non-zero value of Eq. In fact, this leads to a qualitative 
conclusion that the hopping process in ZAEB results in a 
non-zero conductivity also for T^OK. Such conclusion 
can be supported by the following additional arguments. 

According to Ambegaokar, Halperin and Langer^ the 
intrinsic transition rate 7^ for an electron hopping from 
a site i with energy Ei to an empty site j with energy 
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Ej (in the simplest case, when kT is small compared to 
| Ej — Ei | , and | Ej — Ei | is of the order of the Debye energy 
or smaller) is well approximated by the "quantum-limit" 
hopping formula^ 

7ij = ^ e -2afl 4 ,-(B J -E,)/*T for E , > ^ (7) 

7y=7oe" 2aH « forEjKEk, (8) 

where 70 is a constant as defined elsewhere^* and Rij is 
the distance between the centers i and j. To estimate 
the time averaged hopping probability between the cen- 
ters i and j, we consider a bidirectional hopping process 
(i — ► j — v i — ► ...). In a classic semiconductor such a 
process always requires thermal activation - because of 
hops to empty states with higher energy, represented by 
Eq. (7). The situation is qualitatively different in the 
semiconductor, where RPs take place. If we introduce a 
parameter i>ij denoting the probability that due to the 
RP (not due to phonon absorption) an occupied state i 
occurs at higher energy level than an empty state j (i.e. 
Ej > E t before RP, and Ej < E t after RP), then the 
probability of the hopping from site i to site j will be 
Vij7oe~ 2aRij ■ In the case that the ratio 



t lJ j oe ~2aR i:i -{E : ,-E i )/kT ^(Ej-E^/kT V ' 

is greater than 1, the predominating hops at a finite tem- 
perature T are not longer the phonon assisted ones, but 
those "driven" by RPs. At sufficiently low temperatures 
[3ij > 1 for any pair of i and j with nonzero Vij. So that, 
the (bidirectional) intrinsic transition rate at lowest tem- 
peratures (flij >> 1) can be approximated by 

Hj « Hjloe- 2 ^ (10) 

and the probability of electron hops due to dynamical 
changes of the impurity energy levels, ph, in an macro- 
scopic system with RPs can be approximately expressed 
by 

p h <x Prp e- 2aR " . (11) 

Here p rp is an effective parameter describing the prob- 
ability of finding a pair consisting of an occupied state 
and an empty state of less energy at an effective distance 
R* . Such hopping probability represents, in fact, the 
tunneling process across the barrier of "thickness" R*. 
Taking into account a finite value of R* ~ N(Eq)^ 1 / 3 , 
where N(Eq) is given by Eq. (6), the resulting "tunnel- 
ing" conductivity in the ZAEB of non-zero width ^2Eq 
has to be non-zero for T — > K. The crossover from 
the phonon assisted hopping type of conductivity to the 
conductivity driven by RPs is expected at temperatures 
corresponding to /3 ~ 1, where (3 is an effective parameter 
of partial parameters given by Eq. (9). 

The above sketched scenario indicates how the hop- 
ping process in the IB can be induced by "slow" VFs. 
For application of this scenario to SmBg with fast VFs, 



we take into account the following facts. Supposing that 
Fermi level of SmB6 lies in the IB of the width Wjb and 
assuming a subnetwork of impurity states that can be 
considered as ZAEB, the upper estimation for Eq follows 
from the inequality < 2E < Wib < E g , where E g 
is the width of the forbidden gap in SmBg. Such as ex- 
periments give Eg > 2 meV, for purposes of our rough 
estimation we consider Eq to be less than ~1 meV, pro- 
viding ^2 meV for the upper estimation of the energy 
of non-equilibrium excitations. Based on the Heisenberg 
relation we estimate the life time of non-equilibrium ex- 
citations to be longer than the life time of excitations 
with energy 2 meV, t exc = h/2 meV = 3.3 x 1CP 11 s. 
Dynamical changes of the impurity energy levels due to 
Coulomb interaction between the impurity and the sur- 
rounding lattice can be characterized by the time con- 
stant of charge fluctuations. The charge fluctuation rate 
of SmBg estimated from phonon spectroscopy studie a 15 i 16 
is between 200 cm -1 and 650 cm -1 . These values corre- 
spond to the characteristic time between 5.1 x 10~ 14 s and 
1.7 x 10 -13 s, what are the values more than two order of 
magnitude less than the life time of non-equilibrium exci- 
tations. So that the energy of both, occupied and empty 
states can vary within their SIBs for many characteristic 
time constants of charge fluctuations without a change of 
their occupation. This means that SIBs of occupied and 
those of empty states overlap in ZAEB (of the width 2Eq 
centered at Ep). Considering that the dynamics of VFs 
in SmBg is temperature independent down to absolute 
zero, ZAEB persists in the system down to the lowest 
temperatures. Therefore an energy distribution function 
(EDF) of charge carriers at the lowest temperatures has 
to be deviated from the (equilibrium) Fermi-Dirac dis- 
tribution function (FDDF). As a consequence, electrical 
conductivity of SmBg does not converge to zero at lowest 
temperatures, such as only the ground state described by 
FDDF (with all occupied states below Ep and all empty 
states above Ep at T = K) is characterized by zero 
value of electrical conductivity (of hopping type) . More- 
over, taking into account the effect of temperature on 
EDF of charge carriers, the presence of ZAEB of non-zero 
width 2E indicates that the charge carrier subsystem is 
not able to reach temperatures, at which the width of 
thermal broadening, ^kT, is less than 2Eq. This implies 
that the charge carrier subsystem can not be cooled below 
a certain minimum temperature T m i n > 2Eo/k > K. 
(Note that such conclusion also opens a question on the 
relevant definition of the "ground state" in Kondo insu- 
lators.) According to our opinion, just a non-equilibrium 
EDF of charge carriers is a fundamental reason for exper- 
imentaly observed saturation of electrical conductivity of 
SmBg (and other Kondo insulators, e.g. FeSi) 17 i 18 ' 19 at 
the lowest temperatures. 

Based on the discussion above we emphasize that pro- 
posed concept of hopping-type conduction is applicable 
in systems, with energy of localized states driven by lo- 
cal dynamical effects. Moreover, if such effects can be 
induced into the system by external parameters, then 
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these parameters can be used to control the electrical 
conductivity of the system. This indicates a possibility 
to develop new type of active electronic devices. Accord- 
ing to our opinion, a relevant example is the influence of 
local dynamical changes of voltage potential in the gate 
area of a MOSFET on two-dimensional (2D) conductiv- 
ity in the device. An energy difference of two (almost 
identical) localized states in the conducting 2D layer at a 
mutual distance R* is governed by voltage difference be- 
tween the corresponding regions in the gate. Because of 
non-zero electrical resistance between the regions, ther- 
mal noise voltage occurs as follows from Nyquist for- 
mula. This causes a fluctuating energy difference be- 
tween the impurity centers, giving rise to a deviation 
from the equilibrium FDDF and to a hopping process 
(driven by gate-noise voltage), in accordance with equa- 
tion (11). However, in contrast to the above discussed 
valence-fluctuating semiconductors, it has to be taken 
into account that the thermal noise voltage vanishes for 
T — » K. Therefore, if there exists a noise voltage in- 
duced metallic-like behavior [dpjdT > 0) dominating at 
relatively higher temperatures, it has to change to tem- 
perature activated transport (dp/dT < 0) for T — > K, 
and so p — > oo for T — > K. Indeed, studies of cool- 
ing of 2D electrons in silicon MOSFET reported by Prus 
and coworkers^ have shown that some of "metallic" p(T) 
curves turn insulating below ^300 mK3 

Discussion above adverts to a possible role of local volt- 
age fluctuations in metallic behavior in two dimensions^. 
A detail understanding of this phenomenon seems to be 
fundamental, as according to the arguments of Abra- 
hams et al£^ a noninteracting 2D carrier system with any 



amount of disorder has to be localized at zero tempera- 
ture. More light into the problem should introduce ex- 
periments with MOSFET-type devices devoted to careful 
studies of the influence of several types of local dynam- 
ical changes of voltage potential in the gate area on (i) 
the electrical conductivity, on (ii) the noise current, and 
on (hi) the charge carrier temperature in 2D conducting 
layer. Especially the latter point is of special importance 
as it can bring direct experimental proof of the model 
presented in this paper. It would be also reasonable 
to verify, whether a novel design of the gate electrode 
aimed to the controlled generation of local (say nano- 
range) voltage fluctuations or oscillations (for instance 
due to the current flow in the area of nanostructured 
gate) can be considered as another, or even alternative, 
parameter controlling the 2D conductivity in the novel 
MOSFET-type devices. 

In summary, we have provided arguments that dynam- 
ical changes of energy of localized states due to valence 
fluctuations can be the clue for understanding experimen- 
tally observed unusual metallic like transport in Kondo 
insulator SmBg and related systems. The proposed hop- 
ping scenario also introduces new light into the problem 
of anomalous metallic-like 2D conductivity and indicates 
a favorable route for research that can bring experimental 
verification of the discussed phenomena. 
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